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Abstract 

We propose an explanation for the LSND anomaly based on quantum decoherence, postulating 
an exponential behavior for the decoherence parameters as a function of the neutrino energy. 
Within this ansatz decoherence effects are suppressed for neutrino energies above 200 MeV as 
well as around and below few MeV, restricting deviations from standard three-flavour oscillations 
only to the LSND energy range of 20-50 MeV. The scenario is consistent with the global data 
on neutrino oscillations, alleviates the tension between LSND and KARMEN, and predicts a null- 
result for MiniBooNE. No sterile neutrinos are introduced, conflict with cosmology is avoided, and 
no tension between short-baseline appearance and disappearance data arises. The proposal can be 
tested at planned reactor experiments with baselines of around 50 km, such as JUNO or RENO-50. 
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I. INTRODUCTION 


The three neutrino mass and mixing scheme has been established as the standard 
paradigm to explain solar, atmospheric, long baseline and reactor experiment neutrino 
data. There are however some anomalies that cannot be explained within this standard 
paradigm. Arguably the most signihcant one is the LSND anomaly [1]. The canonical 
solution to the LSND anomaly is the four neutrino mixing scheme that includes a sterile 
neutrino with mass of order of 1 eV and a small mixing with v^. and This solution suffers 
from some drawbacks. Most importantly, within this scheme there is a tension between ap¬ 
pearance and disappearance experiments, see Eli for recent analyses. Moreover, a sterile 
neutrino with mass and mixing parameters to solve the LSND anomaly is in tension with 
cosmology Hi- 

In view of these tensions, quantum decoherence has been suggested in the literature to 
explain LSND [UHH]. It is hypothesized that the evolution of quantum states receives a cor¬ 
rection relative to the prediction of standard quantum mechanics. Such effects could arise 
in certain quantum gravity scenarios [9]. To explain the results of short baseline neutrino 
experiments, a phenomenological approach is taken to determine the form and magnitude 
of decoherence. Neutrino oscillation is a quantum interference effect over macroscopic dis¬ 
tances, which can be sensitive to small corrections to quantum mechanics. The idea of 
involving quantum decoherence to explain the LSND anomaly was first proposed in 
As discussed in [8], within the framework proposed in the decoherence effects exceed 

the upper bound from the NuTeV experiment [10] in which neutrinos have an average energy 
of 75 GeV. 

Considering this observation, the so-called soft decoherence scenario was suggested in |S] 
as a solution to the LSND anomaly. Within this scenario, the decoherence effects rapidly 
decrease with neutrino energy, avoiding the NuTeV bound and leaving neutrino oscillations 
in experiments with GeV scale neutrino energies unaffected. Furthermore, decoherence is 
restricted to the 1-3 sector, while the 1-2 sector is not modihed. This explanation does not 
suffer from the appearance-disappearance tension of sterile neutrino models. On the other 
hand, with the power law energy dependence that was assumed in [8], reactor neutrinos 
undergo quantum decoherence after propagating distances of few cm. This implies that no 
oscillation due to Amh would be possible along the distance between near and far detectors 
of reactor experiments. For this reason, the soft decoherence scenario of [8] is now excluded 
by the results of the Daya Bay mm and RENO [T^ IH] experiments]^ In the present 
paper, we revisit the decoherence scenario by modifying the power law assumed in [8] to an 
exponential energy dependence of the decoherence parameter, leading to an explanation of 
LSND consistent with all existing data. 

The outline of the paper is as follows. In Sec. |^we review the decoherence scenario and 
introduce the ansatz for the exponential energy dependence of the decoherence coefficients. 


Sec. Ill contains the numerical results of our analysis of the relevant oscillation data, showing 


that the scenario can explain LSND without being in conflict with other data. In Sec. IV 


^ Note that the DoubleChooz experiment 
comparison is available to date. 


does not (yet) exclude this scenario, since no data on near-far 
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we discuss further implications of the scenario and predictions for future experiments. In 
particular, we show that planned intermediate baseline (50 km) reactor experiments can 
provide a crucial test of the framework. We conclude in Sec. |V| 


II. QUANTUM DECOHERENCE AND THE LSND ANOMALY 


In the quantum decoherence framework, the evolution of the density matrix for neutrinos 
can be described as 

| = ( 1 ) 

where H is the Hamiltonian and T>[p] parameterizes the decoherence effects. Maintaining 
complete positivity leads to the Lindblad form for Dip] (THl II7| 

®[p1 = E [{'’• - 2D„,pDl] (2) 

m 

where Dm are general complex matrices. Unitarity then requires Dm to be Hermitian. If 
we further impose conservation of average energy, we hnd [H, Dm\ =0. As a result, in the 
neutrino mass basis. Dm and H can be simultaneously diagonalized 


H = Diag[hi, / 12 , hs ], Dm = Diag[dm,i, dm, 2 , dm, 3 ], 


( 3 ) 


where hi = + (adopting the equal momentum approximation for the mass states), 

and dm,i are unknown energy dependent real quantities with dimension of Solving 

Eq. ([^, we hnd 


pit) 


Pii(O) pi2(0)e-('^i2-iAi2)t 

P2i(0)e-(^2i-iA2i)i ^2^(0) p23(0)e-(^23-*A23)t 

.P3l(0)e-(^31-iA3l)i ^^^^Q)g-(^32-iA33)t ^33(0) 


( 4 ) 


in which 


hij — ^ ^y dm,i dm,j) and Ajj — hj hi ~ 


Am?, 




2E. 


( 5 ) 


Obviously, jij = jji and Aij = — A^j. This means 'jij is symmetric under hipping i ^ j. In 
the following, we assume that only one term contributes in the sum and we drop the index 
m. The havor conversion probability can be written as 

ij 


where Uai are the elements of the PMNS matrix [18]. The density matrix p'y\t) is given by 


Eq. (4) and pp(0) = PiJ^(O) = UaiU*j. The havor conversion probability for antineutrinos, 
will be given by a similar formula, replacing U with U*. 
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In the soft decoherence scenario of [ 8 ], a power law energy dependence of the decoherence 
coefficients has been assumed, di oc E~'^ (r > 2), suppressing decoherence effects for E > 
100 MeV. However, as mentioned in the introduction, this leads to strong decoherence 
effects at low energies and is by now excluded by Daya Bay and RENO results. In this 
work, we therefore propose a modified energy dependence of the decoherence parameters 
and we conjecture an exponential dependence on energy for di as follows: 


di = exp 



(7) 


where 70 is a constant parameter with dimension of mass, universal for all mass eigenstates. 
Ei are also constant parameters with dimension of mass but can in principle take different 
values for different mass eigenstates. The power n can take any arbitrary number. In 
line with the idea of soft decoherence, we take a value for n and Ei for which at energies 
> few X 100 MeV, the decoherence parameters become suppressed rapidly enough not to 
have any effects at experiments such as MINOS ra. T2K [20], atmospheric neutrinos [2T] 
and etc. In the same way this predicts null-results for short-baseline experiments with 
E > 200 MeV such as MiniBooNE [22[|23] , CDHS [21|, NOMAD [20], NuTeV [10] and etc. 
We found that with Ei < 100 MeV and n > 2, this requirement is fulfilled. Unless it is 
stated otherwise, we take n = 2 for definiteness throughout this paper. 

To avoid constraints from the long-baseline KamLAND reactor experiment [26], we re¬ 
strict the scenario to di ~ ^2 or equivalently to 712 ~ 0 [ 8 ]. In the limit |Ei — E 2 I Ei 
with taking n = 2 and E < Ei, we find 712 ^ 470 exp{-2E‘^ / Ef)E^{E 2 - Eif/Ef. To avoid 
bounds from KamLAND, 712 should be much smaller than ~ (200 km)“^ at E ~ few MeV 
which for Ei 3> MeV means \E 2 — Ei\/Ei (800 km 7 o)“^/^[Ei/(few MeV)]^. At first sight, 
it seems that from solar neutrino data, we can obtain strong bounds on 712 , too. However, 
for long baselines, the interference effects are averaged out and as a result the sensitivity to 
7 jj is lost. This happens for solar neutrinos even before reaching the resonance region inside 
the Sun. From a theoretical point of view, it may be natural to assume that di are functions 
of mass: di = f{mi). From mi ~ m 2 7 ^ m 3 , we then expect di ~ ^2 7 ^ d^. 

In the rest of this paper, we shall take 


7 i 2 — 0 and 7 = 713 — 732 — 7o 


/ 

fEy 


r fEY] 

1 exp 

. [esJ . 

— exp 

[- [eJ \ 


( 8 ) 


with n = 2. Notice that the combination in the parenthesis is less than or equal to 1 and 
hence, 7 < 70 . For E 3> Ei,E 3 , we have di,d 3 —)■ 0 and 7 will therefore exponentially 
converge to zero. For E Ei, E^, 7 will also be small and suppressed by [F'”'(Fj“"' — 

Only for E ~ Ei,Es, the value of 7 can be sizable and decoherence effects can be significant. 
Note that the suppression of decoherence at low energies works only for a universal coefficient 
7o. Hence, the assumption that 70 is independent of the neutrino mass is crucial for our 
scenario. 

For A 21 T <C 1, we can write 


P^e-(7,^) =7^/.e(7,^) = 7"ei.(7,^) -2|t/^3nf/e3p [l - cos(A3iL)] (9) 

Peeil,L) = Peeil,L) ~ 1 - 2\Ue3\\l - \Ue3\^) [l - e-^Sos(A31L)] , (10) 

F^^( 7 , L) = P^^( 7 , P) ^ 1 - 2|P^3p(l - [1 - cos(A3iP)] . (11) 
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For 'jL —)■ 0 the quantum decoherence is turned off and the flavor conversion probability 
becomes equal to that in the standard three neutrino oscillation scenario. 


III. ANALYSIS OF SHORT BASELINE AND REACTOR NEUTRINO DATA 

In this section, we present the results from a numerical analysis of relevant data and 
determine the allowed range of parameters which can account for the LSND anomaly without 
being in conflict with any other experimental results. 


A. Description of the used data and analysis details 

In our analysis, we focus on the LSND electron antineutrino excess events in the energy 
range from 20 MeV to 60 MeV [Tj. We extract the data points as well as the background 
from Fig. 24 of pQ. The data sample shown in that hgure was obtained by applying the 
analysis cut > 10, see section VII-C of [T] for the dehnition of the R-y variable. To predict 
the number of events in each bin within the decoherence scenario, we normalize the total 
number of events for —)• Pg) = 1 to 33300 as indicated in table VIII of [Ij, multiplied 
by 0.39 which is the efficiency of the R-y > 10 cut (see table IX of m)- The is dehned as 
the sum of squares of the difference between prediction (signal+background) and observed 
number of events per bin divided by the square of the uncertainty. The sum is over bins 
with 20 MeV < E < 60 MeV (10 bins) and the uncertainties in each bin are obtained from 
the error bars on the data points in Fig. 24 of [1], which account for both systematic and 
statistical uncertainties. We have checked that our analysis reproduces the allowed region 
for standard oscillations obtained in jl] with good accuracy. 

We also take into account the results of the KARMEN experiment [27], which observes 
15 events in the energy range from 16 MeV to 52 MeV with a predicted background of 
15.8 ± 0.5 events. Any explanation of the LSND anomaly has to address the null-result of 
KARMEN, taking into account the very similar experimental conhguration, with the main 
difference being the somewhat shorter baseline of KARMEN. Again we perform a fit to 
the binned energy spectrum (9 bins) and we can reproduce the official results in terms of 
sterile neutrino oscillations to good accuracy. For short baseline experiments such as LSND 
and KARMEN, A 31 L 1, so we can use Eq. (|^ to write the conversion probability for 
neutrinos and antineutrinos as follows: 

Fpg(7,L) = P,eh:L) = 2\Uy,s\^\U,s\^ (1 - e-^^) ^ (1 - e-^^) . (12) 

As discussed before, the data from Daya Bay and RENO, being consistent with the 
standard three neutrino oscillation scheme, can put bounds on the decoherence parameters 
El and E^. To derive the bounds, we analyze the energy spectrum of the Pe flux at the near 
and far detectors of Daya Bay shown in Fig. 2 of m- We read the data points for near 
detectors (EHl and EII2) and far detector (EH3) from the upper panel in pairs of panels 
shown in Fig. 2 of m, 75 data points in total. We read the background for each detector 
from the inset panels in this hgure. Finally, having extracted the data and background, to 
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calculate the number of events per bin without oscillation {i.e., for P{i'e —t = 1), we use 
the data points displayed in the lower panels [(data-background)/predictions] of Fig. 2 of 
m- To calculate the number of events within the decoherence scenario, we then multiply 
this number with the probability in Eq. (10), averaged over cross section, flux, and energy 
resolution. To compute the we equate uncertainties for each bin to the root of number 
of events per bin (he., the statistical uncertainty). The overall flux normalization is taken 
to be a free parameter to be hxed by the combined near and far detector £t. The distances 
between the various reactor and detector sites of the Daya Bay experiment are taken from 
table 2 of [T2] . 

We also include the spectrum of fifty thousand inverse beta-decay candidate events of the 
far detector of the RENO experiment [13] and compare it with the prediction. Data points 
are taken from the right panel of Fig. 4 of [TT] (26 data points), where the background 
is already subtracted. To compute the prediction of the decoherence scenario for each 
bin we multiply the oscillation prediction shown in Fig. 4 of [T3| by the averaged survival 


probability in Eq. (10) and divide by the oscillation probability with sin^ 26*13 = 0.094 and 
|Am|i| = 2.32 X 10“^ eV^ as stated in Fig. 4 of [14] . These mass and mixing parameters are 
the best fit values that Ref. |E1 derives by using a MC simulation to fit both near and far 
detector data. For values of Ei and E 3 of interest for solving the LSND anomaly {Ei E 3 ^ 
few 10 MeV), decoherence at the near detector is negligible (7L 1) and Eq. (10) converges 

to the standard oscillation formula. Hence, using the far detector prediction based on the 
near detector data (as done for Fig. 4 of [H]) should be a good approximation. Notice, 
however, that including the RENO results does not much change the overall results for the 
decoherence £t, which is dominated by Daya Bay data. 


B. Results of our fit 

Remember that at LSND should be of order fewx 10“^ to account for the observed 
excess. For |Pe3p — 0.02, the value of 7L for LSND should be of order of 0.1 to explain 
the anomaly. From 7 < 70 and L = 30 m, we hnd that 70 has to be of order of 0.01 m“^ 
or larger. Larger values of 70 require Pi ~ P3 to cause partial cancelation, see Eq. (|^. To 
explain the LSND anomaly, we demand that 7 ~ 70 ~ 0.01 m“^ at P ~ 30 MeV and to 
avoid the bounds from reactor experiment as well as from higher energy experiments, we 
require 7 70 for both P 3> 30 MeV and P ~ few MeV. That means Pi and P3 should 

be of order of 10 MeV. Fig. shows 7L versus energy taking typical values for decoherence 
parameters. The left panel of that figure corresponds to a parameter choice close to the best 
fit value of our model. As seen from Fig. 0(1 eft), at 7o = 0.01 m the effect of decoherence 
is negligible for energies above 200 MeV. Thus, the bounds from short-baseline experiments 
such as NOMAD, CDHS, or NuTeV are satished. In other words, like in the soft decoherence 
scenario, the tension between appearance and disappearance experiments plaguing the 3-1-1 
sterile oscillations is solved. Furthermore, the standard oscillation results for experiments 
with 0(1 GeV) neutrinos such as MINOS, T2K, or atmospheric neutrinos are not affected. 
For Pj > 200 MeV (see right panel of Fig. [^, decoherence effects can potentially show up 
in the low energy bins of T2K as well as in the sub-GeV atmospheric neutrino data. 
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FIG. 1. Dependence of jL on energy for different baselines corresponding approximately to LSND, 
MiniBooNE, medium baseline reactor, long baseline accelerator and atmospheric neutrino experiments. We 
have taken n = 2, 70 = 0.01 m“^ for both panels, and Ei = E 2 = 20 MeV, E 3 = 55 MeV {Ei = E 2 = 
60 MeV, E 3 = 200 MeV) for the left (right) panel. 



FIG. 2. Decoherence prediction for LSND for 70 = 0.01 m ^, Ei = E 2 = 18 MeV and F3 = 63 MeV 
compared with data. 


Our main focus is on a range of parameters for which reactor and T2K experiments are 
unaffected. As a result, a reanalysis of Daya Bay and T2K will approximately yield the 
same value for 613 as in the standard oscillation case. We fix the values of the standard 
neutrino parameters (including 6 ^ 13 ) to the best fit value of the global analysis from [2H]- We 
find that within this scenario with n = 2, LSND data can be explained with a satisfactory 
p-value of 68 % with three unknown parameters fitted to 70 = 0.01 m“^, Ei = E 2 = 18 MeV 
and E 3 = 63 MeV. The spectrum of events at LSND for these values is shown in Fig. 
The figure demonstrates that data and prediction of the decoherence scenario are in good 
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FIG. 3. Constrains on the parameters from short baseline and reactor experiments at 90% C.L. 
taking n = 2 and 70 = 0.01 m“^. The region below and to the left of the dotted curves is allowed by 
KARMEN, the region between the dark-blue solid curves is allowed by Daya Bay and RENO, the thin 
magenta curves delimit the regions allowed by LSND. The regions consistent with all data are inside the 
thick red curves, with the cross indicating the best fit point. 


agreement. In the following analysis we will fix 70 to the LSND best £t value of 0.01 m 

Fig. shows the constraints from short baseline and reactor neutrino experiments on Ei 
and E 3 at 90 % C.L., fixing n = 2 and 70 = 0.01 m“^. As expected the bounds are symmetric 
under Ei ^ E^. Due to the exponential dependence on the distance L, the difference in 
the baselines for LSND and KARMEN (30 m versus 18 m) leads to a better consistency of 
the two results than in the case of oscillations. Fig. shows that at 90% C.L. KARMEN 
only marginally constrains the LSND allowed region, compare thin magenta (LSND) and 
dotted black (KARMEN) curves. For Ei ^ 8 MeV, the bound from the reactor neutrino 
experiments becomes stringent and practically only the narrow region with Ei ~ E^ is 
allowed. But for Ei^E^ > 15 MeV, the bounds from reactor experiments are relaxed. The 
thick red curves in the plot show the globally allowed region. The best fit is marked by a 
cross in the plot and it is located at Ei = 20 MeV and E^ = 56 MeV. We have clipped the 
figure at ^ 1,^3 = 200 MeV because for larger values T2K and atmospheric neutrinos will 
also be affected, as visible from the right panel of Fig. 

Table [^shows Xmin degrees of freedom and goodness of fit (GOF or p-value) for various 
short baseline and reactor neutrino experiments. Notice that the p-value for LSND given in 
the table (he., for the case that 70 is hxed and Ei and E 3 are treated as free parameters) is 
better than the aforementioned p-value that we obtain when we treat 70 as a free parameter 
along with Ei and E^. This reflects the fact that for 70 > 0.01 m“^, the minimum value 
of over the Ei and E^ plane does not change much by varying the value of 79 . Let us 
comment on the somewhat large p-value of 93% for the reactor analysis. If Daya Bay and 
RENO are analyzed separately we hnd x^j^/DOF values of 55/72 and 23/24, respectively. 
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Data xLn/DOF GOF xIg/DOF PG 


LSND 

4.8/8 

77% 



KARMEN 

7.0/7 

43% 



Daya Bay and RENO 

78/98 

93% 



LSND-bKARMEN 

14/17 

66% 

2.3/2 

32% 

LSND-bKARMEN-bReactor 

93/118 

96% 

3.2/4 

52% 


TABLE I. 


Xmin/DOF and goodness of fit (GOF) for different combinations of short baseline and reactor 


neutrino data. The last two columns quantify the consistency of different experiments, see Eq. (13) for 
the definition. Ei = E 2 and E 3 are taken as free parameters to fit the data and the rest are fixed to 
7o = 0.01 m“^, n = 2 and sin^ 2^13 = 0.085. 


Hence, the too good fit comes from the Daya Bay analysis. This might be related to the 
accnracy of reading data from the plot. Note however, that onr results are based on Ax^ 
values, which are insensitive to the absolute value of the x^- Furthermore, we can reproduce 
the standard 613 result of Daya Bay with good accuracy. 

Consistency of the combination of various experiments is quantified by the so-called 
Parameter Goodness of fit (PG) [29l [30] defined as 


XPG 


Xtot,] 


J2A 


(13) 


where Xtot min global minimum, the sum over i runs over the different experiments, 

and x?min minima of the experiments separately. As seen from the third and fourth 

columns of table [I], the KARMEN and LSND data are in good agreement with each other 
under the decoherence hypothesis, thanks to the exponential dependence of the transition 
probability on L. Furthermore, the short-baseline experiments LSND and KARMEN are 
also in very good agreement with the Daya Bay and RENO reactor experiments. 

From Fig.[T| it is clear that decoherence effects are strongly suppressed for the MiniBooNE 
baseline of around 500 m and neutrino energies above 200 MeV. Thus, our scenario is con¬ 
sistent with a null-result in MiniBooNE, as observed in the energy range E > 475 MeV 
The low-energy event excess between 200 and 475 MeV inn is not explained since the transi¬ 
tion probability is already highly suppressed in that regime. It is necessary to be consistent 
with T2K and atmospheric neutrino data, which requires us to restrict E 13 to values suffi¬ 
ciently low such that not to affect the standard oscillation behavior seen there. Hence the 
low-energy MiniBooNE excess has to find an alternative explanation. 


IV. PREDICTIONS FOR FUTURE EXPERIMENTS AND POSSIBLE EXPERI¬ 
MENTAL TESTS 

First we mention that the so-called reactor |32] and Gallium [33] anomalies cannot be 
explained in the decoherence framework proposed here. At reactor energies and below, the 
decoherence effects are suppressed so we predict neither a reduced reactor neutrino flux at 
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short baselines, nor a reduced neutrino rate in source experiments at Gallium detectors. 
Those anomalies (which are at the level of 3a) should hnd another explanation in the sce¬ 
nario discussed here. Planned experiments at reactors with very short baselines as well as 
radioactive source experiments should lead to null-results. 

Let us now comment on future accelerator-based experiments. Long-baseline oscillation 
experiments such as NOvA |3l] or LBNF [35] use neutrino beams with ^ 1 GeV. As clear 
from Fig. [^we predict no decoherence effects at those energies and hence such experiments 
should obtain results consistent with standard three-flavour oscillations. The nuSTORM 
short baseline neutrino experiment with an average energy of 3 GeV and a baseline of 2 km 
is proposed to test the 3-1-1 oscillation hypothesis [35]. From Fig. we observe again that 
the decoherence effects for this setup are too small so we predict a null signal for such an 
experiment. If nuSTORM hnds no signal for appearance, the 3-1-1 solution will be ruled 
out but the decoherence solution will still survive, while the observation of an appearance 
signal at nuSTORM would exclude the decoherence solution proposed here. The situation 
is similar also for other short baseline neutrino experiment proposals with neutrino energies 
> few X 100 MeV, see e.g., I37|. For the ESS superbeam [38], with a peak energy of 
E ~ 200 MeV some decoherence effects may start to show up if the Ei^^ parameters are 
not too small. We do predict an appearance signal in LSND-like short baseline experiments 
with energies around 30 MeV; see e.g., [39] . 

A crucial test of our scenario might be possible with reactor experiments with baselines of 
around 50 km, such as the JUNO [ID] IH] and RENO-50 [H] projects. Using the formalism 
of Sec. we obtain for the three-flavour Uf. survival probability 


P- 

ee 


1 — sin^ 2012 sin^ _ 1 2Q 

2 2 

^ sin^ 2013 [cos^ 0i2 cos(A3iL) -h sin^ 0i2 cos(A32L)] . 


(14) 


The main goal of those experiments is to observe the small “wiggles” in the energy spectrum 
induced by the term in the second line of Eq. 0. We see that for yL > 1 those features 
will be suppressed due to decoherence. Fig. (left) shows that for baselines of 50 km, 
yL becomes of order one for E^ ^ A MeV, and hence, the fast oscillations in the survival 
probability may become suppressed. 

We perform a numerical study of this effect by simulating the JUNO experiment, using 
information from HU. We normalize the number of events such that for the default exposure 
of 20 kt detector mass x 36 GW reactor power x 6 yr exposure (4320 kt GW yr in total) 
we obtain 10^ events. The energy resolution is assumed to be 3%y/l MeV/E,y. We perform 
a analysis using 350 bins for the energy spectrum. Several systematic uncertainties are 
included as well as the smearing induced by the baseline distribution of 12 relevant reactor 
cores. Further details of our analysis can be found in [42j . 

In Fig. 1^ (left) we compare the expected spectra for standard oscillations (red solid) to 
the decoherence scenario with parameters close to the best ht point (blue dashed). We 
clearly observe that the 0i3-induced modulation of the spectrum (second line in Eq. (14)) 
becomes suppressed in the case of decoherence for neutrino energies above 4 MeV. Thanks 
to the huge number of events this difference is highly signihcant: the Ay^ between those two 
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FIG. 4. Left: Event spectrum at JUNO for an exposure of 4320 kt GW yr. The red solid curve 
corresponds to standard oscillations with Am|^ > 0, whereas the blue dashed curve shows the spectrum 
with decoherence parameters Ei = E 2 = 20 MeV, E^ = 55 MeV, 70 = 0.01 m“^, n = 2. The shaded 
band indicates the statistical error per 0.02 MeV bin. Right: In the shaded regions, JUNO can distinguish 
the decoherence scenario from standard oscillations at more than 3cr (Ay^ = 9). The red curves show the 
90% C.L. allowed regions from the combined analysis of LSND, KARMEN, Daya Bay and RENO, with the 
cross indicating the best fit point (same as in Fig. [^. 


curves is 33, which means that the no-decoherence hypothesis would be excluded at more 
than 5a (1 DOF, i.e., for fixed parameters). The right panel of Fig. shows the regions 
in the {Ei,E^) plane, where JUNO will be able to distinguish the decoherence scenario 
from standard oscillations at more than 3a (Ay^ = 9). We observe that for values of 
Ei,E 3 > 30 MeV, JUNO loses sensitivity, since the decoherence effects will be shifted to 
higher energies and the reactor neutrino spectrum would be very little modified. We note 
however, that for such large values of Ei, decoherence effects may show up in long-baseline 
or atmospheric neutrino experiments. 

Hence, if JUNO does not find any deviation of the energy spectrum from standard os¬ 
cillations, our scenario would be highly constrained. A dedicated investigation of future 
data from reactor, long-baseline, and atmospheric neutrino experiments would be required 
to determine whether an allowed region survives or not. In our analysis we have fixed the 
exponent in Eq. ([^ to n = 2. If JUNO would obain results consistent with standard os¬ 
cillations, one might also test values n > 2. Increasing n, the decoherence for E < Ei,E‘i 
becomes more strongly suppressed which in turn leads to a faster weakening of the sensitiv¬ 
ity of JUNO to the parameter range allowed by LSND. Such investigations are beyond the 
scope of the present work. Note also that we have fixed 70 = 0.01 m“^. Smaller values of 
7o would not allow to fit LSND, as discussed in the previous section. For larger values of 70 
decoherence effects will become larger at JUNO, increasing the sensitivity. 

Finally let us mention that the scenario presented here predicts also large modifications 
for supernova neutrinos, since the energy range of supernova neutrinos (tens of MeV) is 
close to the LSND energy range, where decoherence effects are important. A discussion of 
decoherence effects for supernova neutrinos can be found in [ 8 ] . 
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V. CONCLUSIONS 


We have revisited the idea of quantum decoherence as a solution to the LSND anomaly 
proposed in [8] taking into account the recent results from the Daya Bay and RENO reactor 
experiments. We assume an exponential dependence of the decoherence parameters on 
neutrino energy as shown in Eq. Q. For a suitable choice of parameters the decoherence 
effects can become suppressed for neutrino energies both below and above LSND energies, 
restricting deviations from the standard three-flavour oscillation scenario to the 20-50 MeV 
energy range. In this way neither standard oscillations of MeV neutrinos from the Sun and 
from reactors are modihed, nor the results for neutrinos with energies greater than 200 MeV 
are affected, as relevant for short and long baseline accelerator experiments and atmospheric 
neutrinos. Moreover the scenario becomes free from the famous appearance-disappearance 
tension that plagues the 3-1-1 sterile neutrino solution for the LSND anomaly. We have 
studied the parameter space in which the LSND anomaly can be explained and constraints 
from various reactor and short baseline neutrino experiments can be avoided. Results are 
shown in Fig. We have found a remarkable agreement between KARMEN and LSND data 
within this scenario. The decoherence solution to LSND predicts no effect in MiniBooNE 
and is hence consistent with the MiniBooNE null-result for > 475 MeV. However, one 
should seek another resolution for the low energy excess observed in MiniBooNE between 
200 and 475 MeV as well as for the reactor and Gallium anomalies. 

The scenario predicts results consistent with standard three-flavour oscillations for most 
of the upcoming long and short-baseline neutrino oscillation experiments. However, reactor 
experiments at baselines of around 50 km such as the JUNO or RENO-50 projects will 
provide a crucial test of the scenario for large part of the parameter space. 
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